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Abstract. In this paper, we show that when two systems of differential equations 
admitting a symmetry group are related by a point transformation it is always possible 
to generate invariant schemes, one for each system, that are also related by the same 
transformation. This result is used to easily obtain new invariant schemes of some 
differential equations. 
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1. Introduction 

In modern numerical analysis, the development of geometric integration has become an 
increasingly active field of research. Geometric integration is related to the development 
of new numerical schemes that incorporate additional structures of the differential 
equation which is being discretized in order to reproduce the qualitative features of 
the continuous solution, [1]. 

For differential equations possessing symmetries, a natural thing to ask for, when 
discretizing such equations, is to preserve as much symmetry as possible. There are 
mainly two approaches of dealing with such a problem. The first one consists of defining 
invariant schemes on fixed lattices, [2, 3, 4], and consider only transformations that do 
not act on the lattices. In order to obtain interesting symmetries, transformations 
acting simultaneously on more than one point of the mesh must be considered. A 
second point of view consists of defining invariant numerical schemes over evolutive 
lattices. Such an approach is used when considering groups of transformations acting 
both on the dependent and independent variables. In the latter case, two different 
methods of addressing the issue are found in the literature. The first one is based on 
the application of the moving frame method, [5, 6], while the other uses a method based 
on the infinitesimal generators of symmetries, [7, 8, 9, 10]. In this paper, we shall be 
concerned only with invariant schemes generated using the last approach. 
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The applications of the Lie groups to discrete equations is quite recent compared 
to its continuous conterpart, which goes back to the work of Sophus Lie, [11, 12]. The 
first articles in this new field appeared at the beginning of the 1990-ties, [13, 14]. The 
goal is to develop the applications of the Lie groups to discrete equations into tools that 
will be as powerful as for differential equations. The interested reader can find a review 
of important results and an extensive bibliography in [15]. 

The aim of this work is to systematize a result obtained in a paper of Dorodnitsyn 
and Kozlov, [16]. In the article it shown that the invariant schemes for the one 
dimensional Burgers' equation in the potential form, w t + \w1 = w xx , are related to 
the invariant schemes for the one dimensional linear heat equation, ut = u xx , by the 
point transformation w = — 2\n(u), the same transformation relating the continuous 
equations. We shall see that this is just an example of a more general result. It is well 
known that when two realizations of Lie symmetry algebras of differential equations by 
vector fields can be transformed into each other by a point transformation, then the same 
transformation will map the two differential equations one to the other. In fact, this 
result is also true for the invariants schemes approximating the differential equations. 
In other words, when to differential equations are related by a point transformation, it 
is possible to define invariants schemes for the two equations that are also related by 
the same transformation. This result has some interesting applications. For example, it 
can happen that the discretisation of a differential equation is easier to find if a change 
of coordinates is performed. If that is so, one can then generate an invariant scheme in 
the new variables and find the desired result by performing the inverse transformation. 
A second direct application consists of obtaining solutions of an invariant scheme from 
known solutions of a related invariant scheme. 

The article can be outlined in the following way. First of all, we recall the algorithm 
for generating invariant schemes for a system of differential equations. It is possible 
to generate symmetry-preserving schemes of systems of ordinary differential equations 
(ODEs) as well as systems of partial differential equations (PDEs) with the algorithm. 
Secondly, we show that when two systems of differential equations are related by a 
point transformation, the transformation will also relate the invariant schemes for the 
two systems. Finally, we apply the result to different types of point transformations. 

2. Invariant Difference Schemes 

Let p ) e X = W p , p e N, be the independent variables and u : X — > U C 

M. q , x \— > u(x) = (v}(x), . . . , u q (x)), the dependant variables of the system of differential 
equations 

A(x,u<- n) ) = 0. (1) 

In the equation (1), u^> denotes all the derivatives of u(x) up to order n with respect to 
all the independent variables. Since we want to generate invariant numerical schemes, 
we suppose that (1) is invariant under a group of Lie point symmetries G, of order N, 
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generated by the vector fields 

Vk = J2 Z k,t ( x > u ^ + E u ) d ^ k — 1, . . . ,N. (2) 

i=l a=l 

The set {V fc } forms a basis of a Lie algebra L. 

The discretisation of (1) consists of sampling in the space of independent variables 
X some points x, each labelled by a set of discrete indices, 

The discretisation of the independant variables induces a natural discretisation of the 
dependent ones: 

u m = (u l m ,... 1 u q m ). (4) 

The construction of an invariant numerical scheme consists of determining a consistent 
way of sampling the points x m and defining the evolution of the solution tl<wi SO clS to 
preserve the symmetries of the continuous system of equations. Since the symmetry 
transformations act on independent and dependent variables, we must let the symmetry 
group G act both on the finite difference equations approximating the system of 
differential equations and the lattice on which the approximation is made. 

To adequately model the discrete problem we consider a system of q + p finite 
difference equations 

E k ({(x m+j ,u m+j )} jeJ ) = 0, 1 < k < q + p, where {0} C J C Z p , (5) 

relating the quantities (x, u) at a finite number of points. The set J in (5) is of finite 
order and serves to identify the neighbouring points of (x m ,u m ). By choice, we suppose 
that the first q equations approximate the system of differential equations while the 
last p equations specify the mesh. The p equations for the lattice do not completely 
determine the lattice in general. In fact, they impose restrictions on it. The number 
of points related with each other in (5) depends on the order of the original differential 
equation and the precision we look for. In the continuous limit, we impose that the first 
q equations of (5) go to the system of differential equations while the p other go to the 
identity = 0. Finally, to have an invariant scheme, we request that the system (5) be 
invariant under the group of transformations generated by (2). 

The procedure for generating invariant finite difference equations from a known 
symmetry group is analogous to the continuous case. As for the continuous case, we 
define a prolongation of the group action, but in a different fashion, [15, 7, 8, 9, 10, 
16, 17]. The prolongation of the group action for the discrete problem is realized by 
requiring that the group acts on all points figuring in (5) 

Definition 2.1. Let G be a group of point transformations acting on the space X xU. 
The discrete prolongation of the group action is defined as 

pr g ■ {(x m+j , u m+j )} jeJ = {g ■ (x m+j ,u m+j )} j€j , \/g G G. (6) 
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In fact, the most general definition of the discrete prolongation of the group action 
would require that the transformation acts simultaneously on all the points of the 
discrete space. Since the system of difference equations (5) only involves a finite number 
of discrete points, we can restrict the prolongation of the action to those points. With 
this definition of the discrete prolongation of the group action, it is straightforward 
to derive an expression for the discrete prolongation of an infinitesimal generator of 
tranformations. 

Definition 2.2. Let V k be vector fields, as in (2) defining a basis of L. The discrete 
prolongation of V k is defined as 

Pr V fc := EE + E (7) 

jeJ i=l a=l 

where C m+j = C(x m+j ,u m+j ) and = (p a (x m+j , u m+j ). 

The method for generating a set of fundamental invariants in the discrete case is 
identical to the continuous one, [19, 20]. The only difference is that instead of using the 
continuous prolongation of the group action we use the discrete prolongation. In the 
discrete situation, an invariant involving the points {(x m+ j, u m+ j)}j e jis a quantity that 
satisfies 

/(pr g{{x m+j ,u m+j )} jeJ ) = I{{{x m+j ,u m+j )} jeJ ), V# e G. (8) 

Hence, given a basis of the Lie symmetry algebra L, (2), we look for the quantities 
I satisfying 

pr V k [I ({x m+j , u m+j } jeJ )} =0, k = l,...,N. (9) 

Using the method of characteristics, we obtain a set of elementary invariants I±, . . . , 1^. 
Their number is given by the formula 

H = dim M — rank Z, (10) 

where M is the manifold that G acts on, i.e. M ~ {{x m+ j, u m+ j}j e j}. So dim M = 
p + q x # J, where # J denotes the order of the set J. Z is the N x (p + q) x # J matrix 

z= ; (ii) 

r p JV,l ^7V,2 ,N,q-l,N,q -, / 

formed using the coefficients of the prolonged symmetry generators (7). Since the 
quantities /i,...,/^ form a basis of elementary invariants, any invariant difference 
equation must be written as 

E(I u ...,I,)=0. (12) 

Equations (12) obtained in this manner are said to be strongly invariant and satisfy 
pr V k [E] = identically, [15]. 
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Further invariant equations can be obtained if the rank of the matrix Z is not 
maximal on some manifolds described by equations of the form E({x m+ j, u m+ j}j e j) = 
and satisfy 

pr\ k [E] = 0, k = l,...,N. (13) 

E=0 

Such equations are said to be weakly invariant, [15]. In practice we usually start by 
computing the invariant manifolds since it can facilitate the computation of the set of 
fundamental invariants afterwards. In order for the system of finite difference equations 
(5) to be invariant under the group of symmetries G, it must be formed out of weakly 
or strongly invariant difference equations and so each equation will satisfy (13) on the 
space of solutions. Many different invariant schemes can be formed using the obtained 
invariants. The only requirement is that the continuous limit of the invariant scheme 
gives back the system of differential equations we are discretizing. 

3. Point transformations for invariant schemes 

In this section, we enunciate the main result of this article. 

Theorem 3.1. Let A(x, u^) = and A(x, u^) = be two systems of differential 
equations related by an invertible point transformation 

ip: X x U ^ X x U, (x, u) h-> (x, u) = (if>\ x (x, u),ip\ u (x, it)), (14) 

that also relates their respective symmetry groups G and G. Then under the point 
transformation (14) invariants schemes of A(x, u^) = are mapped to invariant 
schemes of A(x, u^) = 0. 

Proof. Let 

{h,-..,Q, (15) 

be a basis of discrete invariants under G and 

E k ({(x m+j ,u m+j )} jeJ ) = 0, k = 1, . . . ,p + q, (16) 

a invariant scheme of A(x,u^) = 0. We want to show that under the point 
transformation (14), the invariant scheme (16) is mapped to an invariant scheme of 
the system of differential equations 

A(x,u^) = A(ip\ x (x,u), (iP\ u (x, «))(">) = 0. 

First of all, under the point transformation (14) the set of points {(x m+ j, u m+ j)}j e j 
is mapped to 

{(Xm+j, Um+j)}j€J = {i'iXm+j, U m+ j)}j(zj. (17) 

Hence, a set of elementary invariants of G relating the discrete points (17) is given by 

Ik{{%m+j, u m+j}jej) = Ik({^> 1 (Xm+j,U m+ j)}j € j), k — 1, . . . fX, (18) 
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where the Ik are given in (15). Indeed, using the fact that the relation between the 
symmetry groups G and G is givent by 

9 = ^ogo^-\ VgeG. (19) 

we see that the invariance condition (8) is satisfied by the quantities Ik, 

4(pr g{(x m+j ,u m+j )} jeJ ) = 4({V ; ~ 1 o^ogoip- 1 o ip(x m+j ,u m+j )} jeJ ) (20) 

= Ik\\ x m+ji u m+j}jej) = Ik({ x m+j, u m+j}j£j), V (? G G. 

So all strongly invariant equations of (16) are mapped by (14) to new strongly invariant 
equations in terms of the discrete variables (17). 

The same affirmation is true for weakly invariant equations. Given a weakly 
invariant equation E({(x mj , u m+ j)}j e j) = 0, the equation 

E({(x mj ,u m+j )} jeJ ) = S({V'~ 1 (x mj .,S m+J )}j 6 j) = (21) 

is weakly invariant under G. Indeed, 

E(prg{(x mj ,u m+j )} jeJ ) 

E({(x mj ,u m+j )} je j)=0 

= £({^ _1 o^jogo^j- 1 o^(x m .,W m+ j)W) =0. (22) 

So from (20) and (22) we conclude that the new system of finite difference equations 

E k ({(x mj ,u m+j )} jeJ ) = EkH^ixm^Um+j)}^) = 0, k = l,...,p + q, (23) 

is invariant under the group of symmetry G. 

The convergence of (23) to A(x,u) = 0, in the continuous limit, is verified since 
the diagram 

A(ar, w (n) ) = — A(x, u^) = 



Continuous 



limit Continuous 



limit 



(24) 



E(\x m+ j, u m+ j}j e j) — < E{\x m+ j,u m+ j)}j e j) 



commutes. 



□ 



Also, let us mention that any exact solution u = f(x) of (5), which means that 
E({(x m+ j : f(x m+ j)}j eJ ) = is satisfied identically, is mapped by (14) to an exact 
solution of (23) u = tp\ u o since 

E({(x m+j , u m+j )} jeJ ) = E({ip o i/j-^xm+j, f{x m+j ))} jeJ ) = 0. (25) 
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4. Applications 

In our applications, we restrict ourself to scalar differential equations involving at most 
two independant variables. 

To simplify the writing we introduce the following notation 

( X m,ni X m,ni u m,n) = (x , X ( X m,n±l i X m,n±l i u m,n±l) = ±i % ±i u ±) i (26) 

( x m+l,m x m+l,m u m+l,n) = (x , X , u) , { X m—l,w X m-l,m u m—l,n) = (x , X , u) , (27) 

and introduce the steps 

Ax*± = ±(x i ± - x), Aa* = & - x\ Ax* =x i -x i , i = 1, 2. (28) 

There are a lot of interesting point transformations that can be considered. We 
have chosen to look at three different transformations. The first application is concerned 
with the hodograph transformation. The transformation is used to generate symmetry- 
preserving schemes of new equations and obtain exact solutions from known ones. In the 
second example we consider the wave equation in one spatial dimension with a source 
term. By reformulating the problem in the characteristic variables, we shall see that it 
is easy to derive an invariant scheme. Then, by taking the inverse transformation, an 
invariant scheme in the original system of coordinates is obtained. Finally, we investigate 
a particular example involving a change of variables from cartesian to polar coordinates. 



4-1. The hodograph transformation 

Let us start by recalling the definition of a hodograph transformation for a scalar 
differential equation. 

Definition 4.1. Let p ) G MP, and u(x) : MP — > R. The transformation 

H : R p+1 -> R p+1 

1 2 1 ( 29 ) 

(x , . . . , x p , u) i— > (y = u, x , . . . , x p , v = x ) 

is called a pure hodograph transformation. 

The quantity v now plays the role of the independant variable. Under such a 
transformation, the derivatives transform as 

U x i = , I — 2, ... ,71, 

Vy 

1 ^yy ^yyy o ^ yy 

^x 1 i ^x 1 x 1 / Ux 1 x 1 x 1 J ^ rT 1 

Vy vl v 4 y v 5 y 

and so on, [18]. 

This transformation is found in the study of nonlinear differential equations. Such 
a transformation is usually used to linearize differential equations. In our examples, we 
will be interested in going the opposite direction. Given a linear differential equation 
and its invariant scheme we use what we have seen in Section 3 to obtain an invariant 
scheme for the nonlinear differential equation related to the linear one by (29). 
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4-1.1. First Order Inhomogeneous Linear Equation Let us consider the first order 
linear inhomogeneous ODE 

u x - A'(x)u - B'(x)e A(x) = 0. (30) 
This differential equation admits the two dimensional symmetry algebra 

Vl = e A ^d u , V 2 = (u - B(x)e A ^)d u (31) 
and its general solution is 

u(x) = (B(x) + c)e A{x \ (32) 

where c is an integration constant, [10]. 

To discretize (30) we need one discrete variable, m e Z. The schemes that we will 
generate will involve the minimum number of points necessary to approximate the first 
order derivative. This means that it will only involve two points. 

In order to generate an invariant scheme of (30) we first start by finding a set of 
fundamental invariants on the two-point scheme {(x,u), (x + ,u + )}. Hence we look for 
the quantities that satisfy 

( e A(*) dv + e A ^d u+ )I(x, x +) u, u + ) = 0, 

((« - B{x)e A ^)d u + (u + - B(x + )e A ^)d u+ )I(x,x +: u,u + ) = 0. 

The solutions of (33) are 

h = x, I 2 = x+. (34) 

With these two invariants it is not possible to obtain a symmetry-preserving scheme of 
(30). However, the symmetry algebra (31) admits the invariant manifold 

u + e' A(x+) - ue~ A{x) - B(x+) + B{x) = 0. (35) 

Hence an invariant scheme of (30) is given by the system of two equations 

u + e~ A{ - x+) - ue~ A{x) - B(x+) + B(x) = 0, 

(36) 

h - h = e, 

where e is a parameter that goes to zero in the continuous limit. This system of invariant 
difference equations forms an invariant scheme of (30), since in the continuous limit (36) 
goes to (30). We recall that there are no recipes to obtain (36). The only requirements 
are that the system (36) must be formed out of the elementary invariants (34) and the 
weakly invariant equation (35) and that in the continuous limit we recover the ODE 
(30). 

By substituing the continuous solution (32) into (36) it is immediate to verify that 
it is also an exact solution of the discrete problem [10]. 

Now, if we apply the hodograph transformation H : IR 2 — > M. 2 , (x, u) 1— > (y = u, v = 
x), equation (30) transforms to the nonlinear ODE 

v y (A'(v)y + B'(v)e A{v) ) -1=0 (37) 
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The symmetry generators of (37) are 

Vi = e A ^d y) V 2 = (y - B(v)e A ^)d y . (38) 

If we apply the hodograph transformation to (36) we obtain 

E x = y+e- A{v+) - ye- A(v) - B{v + ) + B(v) = 0, (39a) 
E 2 = v + -v = e, (39b) 

and this system of equations satisfies the infinitesimal invariance condition 

prV fc [£j_ _ =0, i = 1,2, k = 1,2. (40) 

E 1= 0,E2=0 

When taking the continuous limit of (39a) and (396) we recover from (39a) the nonlinear 
ODE (37) while (396) goes to the identity = 0. The invariant scheme obtained is quite 
different from standard discretisation. Indeed, in standard numerical methods, the 
lattice can be chosen to have a variable step size, but the choice does not incorporate 
any information on the problem being discretized. Usually, all the information is 
incorporated in the finite difference equation approximating the ODE. In the invariant 
scheme, (39a) and (396), we have the opposite situation. Solving (39a) and (396) we 
find that 

V = me + V(j an d y = (B(v) + c)e A{v \ (41) 

where e, t>o and c are constants. The second equation of (41) defines the evolution of the 
mesh in such a way that the difference in v between two iterations is constant. Hence, all 
the information of the continuous problem is now incorporated entirely in the definition 
of the mesh. In Fig. 1, we have plotted a particular solution to illustrate the situation. 

4-1.2. One Dimensional Linear Heat Equation The one dimensional linear heat 
equation 

u t = u xx (42) 
admits an infinite dimensional group of Lie point symmetries generated by ,[19], 
Vi = d x , V 2 = <9 t , \ 3 = ud u , V 4 = xd x + 2td t , 

(43 a) 

V 5 = 2td x - xud u , V 6 = Atxd x + At 2 d t - (x 2 + 2t)ud u , 

\ a = a(x,t)d u where a t = a xx . (436) 

To discretize (42), we need two discrete indices, (m,n) G Z 2 . 

To perform our invariant discretisation procedure of the heat equation we consider 
only the six dimensional Lie algebra (43a). Before computing the invariants of (43a), 
we point out the important fact that 

At + = (44) 

is an invariant manifold of (43a), which we will take to be one of our equations describing 
the lattice. Hence, the invariants are to be computed on a grid with flat time layers. 
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Figure 1. Solution of v y (l + 4w 3/2 exp^y/v)) = 2y/v, 

exact solution, 

• discrete solution. 



This restriction is important since we want to be able, at any time iteration, to move 
everywhere in the spatial domain. 

The invariants are computed on a scheme involving the points (x,t,u), (x±,t,u±) 
and {x, t : u). This is the minimum of point necessary to approximate the derivatives 
present in the heat equation and generate an explicit scheme. The set of elementary 
invariants is obtained by solving the system of linear partial differential equations 

pr \ k [I(x, t, u, x + , u + , x_, w_, x, t, u)] = 0, k — l,...,6, (45) 

where the quantity I is supposed to depend on the discrete points involved in the scheme 
and the are given in (43a). The solution of (45), by the method of characteristics, 
gives the set of elementary invariants, [21], 



h 



Ax 



+ 



Ax^ 
Axl_ 



Ax 2 + 

~Zt 



M\2 

— exp 

u/ 



-Ax 
2At 



4At Ax H 
Ax + Ax 



+ Ax_ 
2Ax. 



\Ax_ 



ln(^ 

u 



(46) 



At 



Ax + + Ax_ 



u + \ Ax + 
u J Ax_ 



In 



il- 
ii 



The number of elementary invariants is coherent with the formula (10). Indeed, 
since the scheme involves the points (x,t,u), (x±,t,u±) and (x,t,u), the manifold M 
introduce in equation (10) is given by M ~ {x, t, u, x + , u + , x_, tt_, x, t, u} so dim M=10. 
Furthermore, since the prolongation of the 6 vector fields given in (43a) are independent 
we get from (10), \i = 10 — 6 = 4. 
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From the set (46) and the weakly invariant equation (44) we define the explicit 
invariant scheme 



h = 4/ 3 , At + = 0, h = 0, 
which gives in terms of the original variables 



G) 2exp 

At + = 0, 
Ax 



Ax 
2At 



AAt 



Ax, + Ax_ 1 Ax 



In 



V u 



In 



2At 



Ax+ + Ax_ 



In 



(47) 

(48 a) 
(486) 

A./- V -u 7 Ax + 1 ' ' 
Notice that the invariant scheme is not linear even though the PDE is. This is due 
to the fact that we have neglected, in our derivation of the discrete invariants (46), the 
infinite dimensional symmetry generator (436), which states that to any solution of the 
heat equation we can add another solution. 

Let us prove that the system of equations (48a), (486) and (48c) is a valid 
approximation of the heat equation (42) by computing its continuous limit. 

We first start by taking the continuous limit in the discrete variable n, meaning 
that Ax± and At + go to zero. By doing so, the equations (48a), (486) and (48c) go to 

2- 



u 



u 



cxp 



Ax 
2At 



1 - 2At 



V>xx 
U 



= 0, 



u. 



Ax = -2 At— , 



u 



(49a) 

(496) 
(49 c) 



respectively. The development of (49 a) in Taylor series in term of Ax and At gives, if 
we keep only the terms that don't go to zero in the continuous limit, 

la 2 ,m t - /11 \ 2 



-2- 



a u 



r u 



u 



2r 2 



-2- 



u 



+ 2 



u. 



u 



(50) 



Using (49c), we replace the occurances of Ax by At in (50) and take the limit At and 
Ax to zero. Doing so, the equation (50) goes to the heat equation while (49c) goes to 
= 0. 

It can be shown, [21], that 

u(x,t) = — 7==( U . / = T //; + ttj. 



'Ant 



u(x, t) = Ke 



—cx+c 2 t 



t = t m + to, 



X 



X 



(hn + x )(t m + t ), 
h n + xq + 2c(r m + to), 



(51) 



where K, h, Xo, t r and c are constants, are non-trivial exact solutions of the invariant 
scheme composed of the equations (48a), (486) and (48c). 

Now let us perform the pure hodograph transformation H : R 3 — > R 3 , (x,t,u) 1— > 
(y — u, t, v — x). Under such a transformation, the heat equation maps to the nonlinear 
equation 
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The symmetry generators of this equation are, [20] 

Vi = d v , V 2 = d t , V 3 = yd y , V 4 = vd v + 2td t , 
V 5 = 2td v - vyd y , V 6 = 4tvd v + 4t 2 d t - [v 2 + 2t)yd y , 

V a — a(v,t)d y where a t — a vv . 



(53a) 



(536) 



Clearly, since t is unaffected by the hodograph transformation, the expression At + = 
will remain an invariant manifold of the new symmetry algebra. The discrete invariants 
on such an invariant manifold are obtained directly from (46), we just perform the 
hodograph transformation on the set of invariants. Hence we obtain, [21] 



h 



v 



v -v- At 
(v+ — v) 2 (v + — v) 2 



exp 



[v — v) 



2At 



4At 



V-l. — V- 



V + — V 



111 ( — I 



y 



V — V- 



~ (v + — v)(v — v) 2{v + — v) [v — v. 
h = h 



At 



Inl^l 

y 



y 

V 



(54) 



v — V. 



ln(^ 

y 



v + — V 

Using the same invariant expressions as for the discrete heat equation, i.e. 

7 2 = 4J 3 , At + = 0, h = 0, 
we get the explicit invariant scheme for (52), 



y . 

- exp 

y. 



(v — v) 
2At 



At + = 0, 

v — v 2 



At 



v+ 



v+ 



1 



In 



AAt 



v + — V- l v + — V 



1^1 

y 



v — V _ 



ln[^ 

y 



(55) 

(56a) 
(566) 
(56c) 

y J v+ - v 

The computation of the continuous limit of (56a), (566) and (56c) is similar to the 
case of the heat equation. Firstly, by letting the steps generated by the discrete variable 
n go to zero we get respectively for equation (56a), (566) and (56 c) 

2 



V- I V — V- 



y- 



V — V- 



ln 



y± 
y 



y , 

— ex p 



,y, 

= 

V — V 



2At 

y 2 v 2 



= 1 + 2 At ( 

yvy 



1 



y 2 v 2 



(57a) 

(576) 
(57c) 

At yv y 

By developping the two nontrivial equations in Taylor series with respect to the step 
generated by the discrete variable m we get 
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By replacing the appearances of Ay in terms of At in (58 a) with the help of (586) and 
then taking the limit Ay and At to zero we recover from (58a) the PDE (52) while (586) 
goes to = 0. 

Finally, by a direct substitution we can verify that the exact solutions of the linear 
heat equation, (51), are again solutions of the invariant schemes of the nonlinear partial 
differential equation (52), after performing the hodograph transformation on them. 
Namely, we have that 

v — (h n + v )(t m + to), t = r m + t , y = e --» 2 /4^ 

VAnt 

v = hn + v + 2c(r m + t ), t = rm + t , y = Ke~ C0+ ^ >t , (59) 

are exact solutions of the invariant system composed of equations (56a), (566) and (56c), 
where K, h, v , t , r and c are constants. 

From this last example, we notice that even though a pure hodograph 
transformation can linearize a differential equation, it does not mean that the 
transformation on a discrete invariant scheme will have the same effect. Indeed, the 
inverse hodograph transformation can be use to transform the nonlinear equation (52) 
to (42). However, since the invariant scheme of the latter is not linear, the discrete 
invariant scheme for Vt = v~ 2 v yy will not be linearized. 

4-2. (1+1) dimensional wave equation with a source term 

In this section we obtain a symmetry-preserving scheme of the (1+1) dimensional wave 
equation 

u tt -u xx = -4F(«), (60) 

where F{u) is an arbitrary function of u. The symmetry algebra of (60) is generated by 

V 1 = td x + xd t , V 2 = d x , V 3 = d t . (61) 

instead of finding an invariant scheme in the space {x, t, u} we make the change of 
coordinates to the characteristic variables 

^:R 2 ^M 2 , (x,t)^(y,z) = (x + t,x-t). (62) 
Under the transformation (62), the PDE (60) is mapped to 

u yz = F(u). (63) 
A basis of the Lie symmetry algebra of (63) is 

V l =yd y -zd z , ,V 2 = d y , V 3 = d z , (64) 

where V 2 = 1/2(#(V 2 + V 3 )) and V 3 = 1/2(#(V 2 - V 3 )). 
Before computing any invariants we have that 

y-y = 0, z+- z = (65) 

are two weakly invariant equations of the symmetry algebra (64) since they verify the 
condition (13). By including these two equations in our invariant scheme we can compute 
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Figure 2. Invariant lattice of equation (63). 



the invariants on an orthogonal lattice, figure 2. A basis of elementary invariants in the 
space of discrete points {(y, z, u), (y, z, «), (y, z, u), (y+, z, «+), z, u+), i+, u+)} 
is 

Ji = u, I 2 = u + , h = u-, h = u,h = u, I 6 = u+, 
j Ay + j Ay + ~ Az ~ K (66) 

Ay+ 

where Ay + = y + — y. 

From this set of invariants we generate a symmetry-preserving scheme by setting 

~ = F(/i), (67a) 

-'8-' 10 -'lO 

y-y = 0, z + -z = 0. (67b) 

The equations for the lattice, (676), can be solved and give 

y = e(m) n + y (m), z = S(n) m + z (n). 

From (66) we see that it is also possible to impose Ay + = Ay_, by setting I 7 — 1 and 
Az = Az by setting 1 = 1 without loosing any symmetries. So if we do so, the final 
solution for the lattice is 

y = e n + y , z = -5 m - z , (68) 

where e, y , 5 and z are constants. On this rectangular lattice, the equation (67a) 

becomes in term of the discrete variables 

U+ — u — u+ + u , . , . 
~ 5 = F W, (69) 

which correspond to a standard discretisation of (60). The interesting thing to do now 
is to come back to the original variables {x, t, u} by applying the inverse transformation 
of (62) 

r 1 (x,t) = ( 1 -(y + z)^(y-z)). (70) 

First of all, the invariant equations (65) become 

Ax = -At, Ax + = At+. (71) 
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Figure 3. Invariant lattice of equation (60). 



The equations (71) are weakly invariant equations of (61) and are chosen to be part of 
the invariant scheme of (60). The basis of invariants (66) on the invariant manifold (71) 
becomes 

h=u, I 2 = u+, h = u_, h = u, h = u, I 6 = u+, 

Ax + , Aau At . . ( 72 ) 



h = 



Ax + At . /■ 

I 8 = — -±, h = -r:, ho = -4Ax + At, 
Ax + At 



11 ~ Ax.' 

under the transformation (70). By setting the same combination of invariants as in 
(67a) and (67b) and adding the equations I7 — 1 and 1$ = 1 we get an invariant scheme 
of (60) 

h-h J2-/1 ( j 

(73) 



-^8-^10 1 10 
Ax = -At, Ax + = At+, I 7 = 1, I 9 = I. 

Which gives in term of the original variables 

± — - = -4F(u), (74a) 

Ax + At Ax + At 

Ax = -At, At = At, (746) 

Ax + = Ax_, Ax + = At + , (74c) 

where Ax + The solution for the lattice is obtain from the known solution for 

the mesh in y and z and using the transformation (70) 



x = -(e n + y - 5 m - z ), t = -(e n + y + 5 m + z ) 



(75) 



The solution for the mesh (75) is drawn in figure 3. 

We compute the continuous to show that the finite difference equation (74a) is a 
valid approximation of the partial differential equation (60) on the lattice given by (75). 
Taking into account that the lattice implies 

Ax + = At + and Ax+ = Ax + , (76) 
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we find by taking the Taylor series in the steps involving the variation of the discrete 
variable n of (74 a) that 

u x + u t + ^^(uxx + u tt + 2uxt) + 0(Ax 2 + ) 

( Ax Ml (77) 

- \u x + u t + -^{u xx + u tt + 2u xt ) + 0{Ax 2 + )\ | s = -4F(u), 

where the hat over the derivatives of u means that they are evaluated at (x,t). By 
developping the hat derivatives around (x, t) we get 

u tt -u xx + 0(Ax+,At) = -4f(u), (78) 

which goes to the desired limit when Ax + and At go to zero. 

The invariant discretisation (74a) of (60) is quite different from the standard 
discretisation 

u+ — 2u + U- u — 2u + u . « 

f S2— = 

x = e n + xo, t = 5 m + t . 

The difference is due to the counterclockwise rotation of 45 degrees of the rectangular 
lattice in the invariant case. 

4-3. Polar coordinates transformation 

A more interesting transformation of coordinates to be considered is the passage for 
cartesian to polar coordinates. Unlike the other examples the two spaces of variables 
are different since the metric is not the same for each space. 
We consider the elliptic partial differential equation 

u xx + u yy = F(u,x 2 + y 2 ), (79) 

where F is an arbitrary function. The only symmetry of (79) is the invariance under 
rotation 

V = -yd x + xd y . (80) 

Nevertheless it is not obvious how to discretize (79) in cartesian coordinates and preserve 
the rotationnal symmetry. The natural thing to do is to pass to polar coordinates 

x = r cos 9, y = r sin 9. (81) 

In polar coordinates the equation (79) becomes 

u„. + ^u r + -^u ee = F(u, r) (82) 
and (80) becomes the generator of translation in 9 

V = d e . (83) 
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Figure 4. Invariant lattice of equation (79) , 
• discrete points involved in (86), 
O other discrete points. 



Without going through the whole algorithm for finding an invariant scheme, it is 
clear that the standard discretisation of (82) 

u + — 2u + m_ lu + — u 1 u — 2u + u 

e 2 r e r 2 o 2 (84) 

f - r = 0, r + - r = e, 9 + - 9 = 0, 6-6 = 8, 

is invariant under translation in theta since it only involves the difference between two 
values of theta. The solution for the lattice is given by 

r = e n + r , 9 = 5m + 9o. (85) 

where e, 5, ro and 9q are constants. 

Now, if we go back to cartesian coordinates 

r-^W, ^c ta „(!), 

the invariant scheme (84) becomes 
u + — 2u + -u_ 1 u + — u 1 u — 2u + u 



e 



y^Tf e r 2 ■ y 2 5 2 



+ zrrzTi — it — = ( M ' V^+y 2 



^x 2 + y 2 - ^x 2 + y 2 = 0, ^x 2 + + y 2 + - v /x 2 + y 2 = e, (86) 
arctan ^— j — arctan ^— j = 0, arctan (^j — arctan ^— j = 5, 

and is invariant under (80). From (81) and (85) we have that the lattice in the x, y 
variables is given by 

x = (e n + r ) cos(<5 m + 9 ), y = (e n + r ) sin(5 m + 9 ). (87) 

In the Fig. 4 we have illustrated the lattice. By carefully taking continuous limit of (86) 
on the lattice (87) we obtain, as expected, the PDE (79). 
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Conclusion 

It is a well known fact that if the symmetry groups of two systems of differential 
equations are related by a point transformation then these two systems are mapped 
into each other by the same point transformation. In this work, we have shown that 
same result is true for invariant schemes. That is, given a system of differential equations 
and its invariant scheme, any point transformation will map the system of differential 
equations to a new system and the invariant scheme to a new invariant scheme of the 
new system. This result has been used to derive new invariant schemes of differential 
equations from known ones, section 4.1. We have obtained exact solutions of these new 
schemes by applying the point transformation to known solutions. In section 4.2 and 
4.3, we have seen that a change of coordinates can be used to facilitate the computation 
of an invariant scheme of a differential equation. 

The validity of the results obtained in this article rely heavily on the point aspect 
of the transformation ip. The result does not apply to nonlocal transformations 
since such transformations will inevitably break the local character of the symmetry 
t r ansf or mat ions . 
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